Quantitative Determination of Eliashberg Function and Evidence of Strong Electron 
Coupling with Multiple Phonon Modes in Heavily Over-doped (Bi,Pb)2Sr2Cu06+5 
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Super-high resolution laser-based angle-resolved photoemission spectroscopy measurements have 
been carried out on a heavily overdoped (Bi,Pb)2Sr2Cu06+6 (Tc ~ 5 K) superconductor. Taking 
advantage of the high-precision data on the subtle change of the quasi-particle dispersion at different 
temperatures, we develop a general procedure to determine the bare band dispersion and extract the 
bosonic spectral function quantitatively. Our results show unambiguously that the ~70 meV nodal 
kink is due to the electron coupling with the multiple phonon modes, with a large mass enhancement 
factor A ~ 0.42 even in the heavily over-doped regime. 

PACS numbers: 74.25. Jb,71.38.-k,74.72.Gh,79.60.-i 



The interaction of electrons with phonons and other 
collective excitations (bosons) in sohds dictates funda- 
mental physical properties of materials [1]. In a Fermi 
liquid picture, such interaction can be described by the 
Eliashberg spectral function a^F{ui), which accounts for 
the bosonic modes involved in the coupling and their 
strengths. In the conventional superconductors, extrac- 
tion of the Eliashberg function played a key role in pin- 
ning down the phonon as the glue for the electron pairing 
that gives rise to the BCS superconductivity [2]. For high 
temperature cuprate superconductors and other complex 
compounds, one expects that the extraction of the Eliash- 
berg function will also be important in understanding 
the exotic physical properties and the mechanism of high 
temperature superconductivity [3-8]. 

The electron-boson coupling gives rise to the band 
renormalization and the change of the quasiparticle scat- 
tering rate, which can be described by the real part 
(ReS(A:,a;)) and the imaginary part (Iml](fc,a;)) of the 
electron self-energy S(A;,w), respectively. With the dra- 
matic improvement of resolution, angle-resolved photoe- 
mission spectroscopy (ARPES) has emerged as a pow- 
erful tool in probing such many-body effects [5]. Under 
the sudden approximation, it measures the single par- 
ticle spectral function A{k,oj) = (l/7r)ImS(A:, — 
el - Rel](fc,a;)]^ + \hnT,{k , , which is directly re- 
lated to the electron self-energy [5]. ARPES has been 
extensively employed in investigating many-body effects 
in conventional metals [9, 10], graphene [11] and com- 
plex materials [5, 12]. However, a long-standing issue 
with the ARPES technique is the ambiguous "bare band" 
£^[13], which hinders a quantitative determination of the 
electron self-energy ReS(k,e)=e-e^. and the underlying 
bosonic spectral function a^F{uj) [6-8]. 



Coming specifically to high- Tc cuprate superconduc- 
tors, although there is a general consensus on the 
electron-boson coupling as the origin of the universally 
observed 50~80 meV dispersion kink along the (0,0)- 
(7r,7r) nodal direction [14-20], it remains under debate 
on the nature of the boson(s) involved, mainly between 
phonons [15] and magnetic origins [19, 21, 22]. Theoret- 
ical calculations suggest that the electron-phonon cou- 
pling in cuprates is too weak to account for the nodal 
kink [23, 24]. To eventually resolve the issue, it is crucial 
to extract the bosonic spectral function quantitatively so 
that a direct comparison between the experiments and 
the calculations could be possible [6-8]. 

In this paper, we introduce a new approach which can 
unambiguously determine the bare quasi-particle disper- 
sion, and hence enables the quantitative determination 
of the bosonic spectral function. The method utilizes the 
subtle change of the quasi-particle dispersion induced by 
the temperature. This is made possible by the super- 
high resolution laser-based ARPES measurements, which 
present data of the significantly higher quality feasible 
for a quantitative analysis of a modified maximum en- 
tropy method [6]. By carrying out such a procedure 
on a heavily over-doped (Bi,Pb)2Sr2Cu06+5 (Pb-Bi2201, 
Tc~5 K), we are able to determine the bare band and ex- 
tract the bosonic spectral function quantitatively along 
the nodal direction. Our results present a conclusive ev- 
idence that the electron-boson coupling is strong (with a 
mass enhancement factor A « 0.42) even in the heavily 
over-doped cuprate compounds, and it originates from 
the electron coupling to the multiple phonon modes. 
Moreover, our approach is general enough to be appli- 
cable for many other materials. 

The angle-resolved photoemission measurements have 



2 




Momentum k (Tr/a) 



E-Ep(eV) 



k-k,,(A ) 



-0.06 -0.04 -0-02 
E-Ep(i!V) 



FIG. 1: Fine structure in tlie nodal dispersion of a heavily 
over-doped Pb-Bi2201 {Tc - 5 K) measured at 15 K. (a) Pho- 
toemission image taken along the (0,0)-(7r,7r) nodal direction 
with the location of the momentum cut shown in the inset, 
(b) Corresponding photoemission spectra (energy distribution 
curves, EDCs). The EDC at the Fermi momentum kp (red 
curve) shows a dip feature at ~ 73 meV, as indicated by an ar- 
row in the expanded spectrum, (c) Quasi-particle dispersion 
obtained by fitting MDCs (momentum distribution curves) at 
different binding energies. The inset shows the correspond- 
ing MDC width (FWHM) . (d) Energy difference between the 
measured dispersion and a straight line connecting [Ep) and 
—0.3 eV. A peak at ~73 meV is clearly observed as indicated 
by an arrow, (e) Energy difference between the measured dis- 
persion and a straight line connecting and —0.07 eV. Two 
features can be identified at ~ 41 meV and ~ 17 meV. 



been carried out on our newly developed vacuum ultra- 
violet (VUV) laser-based ARPES system [25]. The pho- 
ton energy of the laser is 6.994 eV with a bandwidth of 
0.26 meV. The overall instrumental energy resolution is 
set at 1 meV, which is significantly improved from 10^20 
meV of previous measurements [6-8, 14-20]. The angu- 
lar resolution is ~ 0.3 degree, corresponding to a mo- 
mentum resolution '^0.004 at the photon energy of 
6.994 eV, more than twice improved from 0.009 at 
a usual photon energy of 21.2 eV for the same angular 
resolution. The heavily overdoped (Bi,Pb)2Sr2Cu06+5 
(Pb-Bi2201) single crystals with aTc ^ 5 K were grown 
by the traveling-solvent floating zone method. The sam- 
ples are cleaved in situ in vacuum with a base pressure 
better than 5 x 10"" Torr. 

Figure 1(a) shows the raw data of photoemission im- 
age taken on a heavily overdoped Pb-Bi2201 (Tc~5 K) 
sample along the (0,0)-(7r,7r) nodal direction at 15 K. 
Clear kink is observed near 70 meV in the dispersion 
extracted from MDC (momentum distribution curve) fit- 
ting (Fig. 1(c)), with a corresponding drop in the MDC 
width (inset of Fig. 1(c)). In order to reveal the possible 
fine structures in the dispersion, we plot the energy differ- 
ence between the measured dispersion and a featureless 
straight line connecting (Ep) and —0.3 eV on the dis- 
persion (Fig. 1(d)). A prominent peak shows up near 
^73 meV. In addition, in the energy difference between 
the measured dispersion and another straight line con- 
necting and —0.07 eV (Fig. le), one can also identify 



low energy features at 41 meV and ~ 17 meV, signified 
by the slope changes. We note that these three features 
at ^73meV, ^41meV, and ~17meV are robust and have 
been reproduced in several independent measurements. 
Moreover, a clear dip at 73 meV (Fig. 1(b)) can be ob- 
served on the EDCs (energy distribution curves) near the 
Fermi momentum (kp). 

The underlying bosonic spectral function a^F^uj) can 
be extracted from the measured real part of the electron 
self-energy [6] . The real part of self-energy is related to 
a^F{uj) by: 



ReS(e;T) = 



kT' kTl 



(1) 



where K{y,y') = dx f{x-y)2y' /{x^ -y'^) with f{x) 
being the Fermi distribution function [1] . In the ARPES 
measurements, the real part of the electron self-energy 
is determined from the measured MDC dispersion Cfc by 
Rel](efc) = e/c — Efc, where e° is the bare quasi-particle dis- 
persion for an ad-hoc system with the electron-boson cou- 
pling turned off [26] . Attempts have been made to invert 
the "effective" bosonic spectral function a'^F{uj) from the 
measured MDC dispersion by assuming an empirical bare 
band of a quadratic form e° = a{k — kp^ -\-h{k — kpY [6- 
8]. It was found that the qualitative features of the 
bosonic spectral function do not sensitively depend on 
the choice of the bare band. However, the arbitrariness 
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FIG. 2: Temperature dependence of the nodal dispersion 
measured on a heavily over-doped Pb-Bi2201 (Tc~5 K). The 
dispersions near the kink region are expanded in the up-left 
inset. The bottom-right inset shows MDC width measured at 
25 K, 100 K and 150K. 
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in choosing the appropriate parameters for the bare band 
prevents a quantitative determination of the bosonic 
spectral function [6-8]. Here we propose a general ap- 
proach to circumvent this problem by making use of the 
self-energy change at different temperatures. We assume 
that a^F^ui) and the bare band have negligible tempera- 
ture dependence over a relatively small temperature win- 
dow, and the Eliashberg function can be extracted from 
the self-energy difference at different temperatures by: 



/•OO 

ReS(e; Ti) - ReS(e; T2) = / dw a^Fiuj) 

Jo 



K 



kTi ' fcTi 



K 



(2) 



The successful application of the approach relies crit- 
ically on the high quality of data, as it utilizes a 
small temperature-induced change between two disper- 
sions subjected to noises. Such a stringent requirement 
on the data quality can now be met by the super-high res- 
olution laser-ARPES measurements(Fig. 2). In this pa- 
per, we test our new approach on the heavily overdoped 
Pb-Bi2201, in which the electron-electron correlation is 
relatively weak, and Eq. (2) is believed to be applicable. 

Figure 3 shows the detailed procedure to extract the 
bare band dispersion and the bosonic spectral function. 
To obtain the bosonic spectral function a^F(a;) from 
Eq. (2), one needs to know the difference of the real part 
of self-energy AReE(e) = ReS(e; Ti) - ReS(e; T2) at two 
different temperatures, as schematically shown in Fig. 
3(a). Here, we use two dispersions measured at 25 K and 
100 K. Note that AReE(e) depends on the selection of 
the bare band dispersion. As the result, an iteration ap- 
proach is needed, as detailed in the following: (1) Deter- 
mine AReI](e) based on the current selection of the bare 
band dispersion. The initial value (AReSo(e)) is taken 
as the direct difference between the two dispersions, as 
schematically shown in Fig. 3(a). The result is shown 
in Fig. 3(c); (2) Extract the "bosonic spectral function" 
SF from AReI](e) by Eq. (2), using a maximum entropy 
method (MEM) similar to that in Ref. [6] but with the 
modified integral kernel. The result is shown in Fig. 3(d). 
Note that the significant change of the real part of the 
self-energy is limited to the energy up to ^ 100 meV, 
above which the self-energy difference fluctuates around 
zero and is most likely due to the data noise, as shown 
in Fig. 3(c). Therefore, we assume a maximum bosonic 
mode energy 100 meV when doing the MEM analysis; 
(3) Calculate the real part of self-energy at 25K from the 
extracted bosonic spectral function using Eq. (1). There- 
suit is shown in Fig. 3(d); (4) Determine the bare band 
dispersion by taking the difference between the measured 
dispersion and the calculated real part of the self-energy, 
both at 25K. The result is shown in Fig. 3(b). These 
steps iterate until the result converges. We find that the 
procedure converges quickly and only three iterations are 
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FIG. 3: Extraction procedure of the bare band and the Eliash- 
berg spectral function, (a) Schematic plot for the two disper- 
sions and a bare band to show the real parts of self-energies 
and their difference for a given energy e. (b) Measured nodal 
dispersions of heavily overdoped Pb-Bi2201 at 25K and lOOK 
and the extracted bare bands for different iterations, (c) Self- 
energy difference between the 25K and lOOK dispersions for 
different iterations, (d) Extracted Eliashberg spectral func- 
tion for different iterations, (e) Real part of electron self- 
energy calculated from the extracted spectral function at 25 
K for different iterations. 



needed for this analysis. 

Fig. 4 summarizes the final results of the bare band 
(Fig. 4(a)), Eliashberg spectral function (Fig. 4(b)), and 
electron self-energy (inset of Fig. 4(a)). Although the re- 
sults are obtained by using the two representative disper- 
sions at 25K and lOOK, the self-energy calculated from 
the extracted Eliashberg function (Fig. 4(b)) at other 
temperatures (here we show 60K data) matches very well 
with the measured values, as shown in the inset of Fig. 
4(a), indicating the internal consistency of our approach. 
The obtained bare band is close but less steeper than 
the one given by the LDA calculations [27] (Fig. 4(a)), 
and is significantly different from the empirical one used 
previously (green line in Fig. 4a) [8]. 

Sharp features are clearly identified in the extracted 
Eliashberg spectral function (Fig. 4(b)). Five peaks at 
(14-17), (30-32), (41-44), (54-58) and (70-74) meV can be 
identified below 100 meV, and they are reproducible from 
several independent measurements. In fact, the features 
(14-17), (41-44) and (70-74) meV are clearly visible in the 
raw data (Fig. 1). These sharp features show reason- 
ably good agreement with the phonon modes observed in 
the same material by infrared reflectance (Fig. 4(c)) [28] 
and Raman scattering (Fig. 4(d)) [29]. Note that these 
optical measurements only see phonon modes with zero 
wavevector, and the full band of phonons will exhibit dis- 
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FIG. 4: Extracted bare band dispersion, self-energy and 
Eliashberg spectral function of heavily over-doped Pb-Bi2201. 
(a) Measured nodal dispersion (empty blue circles) at 25 K 
and the extracted bare band (thick red line). For comparison, 
the bare band from the band structure calculations (dashed 
grey line) [27] and the bare band used by previous maximum 
entropy method (thin green line) [27] are also shown. The 
bottom-right inset shows the measured electron self-energy at 
different temperatures and their comparison with those cal- 
culated from the extracted Eliashberg spectral function, (b) 
Extracted Eliashberg functions a^F{uj). (c) Loss function 
from infrared measurements on optimally-doped (red line) 
and over-doped (blue line) Bi-2201 [28]. (d) Raman spectra 
of Bi2201 under different polarization geometries [29]. 



persion in the Brillouin zone, in particular for the higher 
energy phonon near 70~80 meV [30]. Considering the 
rather weak magnetic signal [31], and the absence of mag- 
netic mode at 15 K for the heavily overdoped Pb-Bi2201 
sample (Tc~5 K), our present results leave no doubt that 
the ~ 70 meV dispersion kink is due to electron coupling 
with the multiple phonon modes. 

With the Eliashberg function extracted, it is possible 
to determine quantitatively the strength of the electron- 
boson coupling. We calculate the mass enhancement 
factor X=2j^ ^a^F{uj), and obtain A « 0.42. As 
discussed above, the contribution is mainly from the 
electron-phonon coupling. Our study presents a conclu- 
sive evidence that the electron-phonon coupling is strong 
even in the heavily over-doped cuprate samples. It is 
much stronger than the value suggested by the first prin- 
ciples calculations (0.14^0.20) [23]. Our present work 
thus asks for a re-examination of theoretical calculations 
on the strength of electron-phonon coupling in these com- 
pounds. 

In summary, by taking advantage of the high quality 
data of the laser-ARPES measurements, we have devel- 
oped a new approach to determine the bare band dis- 
persion and extract Eliashberg spectral function quan- 
titatively. The method is general and can be applied 
to many other materials, thus promoting ARPES tech- 



nique to be quantitative in probing many-body effects. 
Our results on the heavily overdoped Pb-Bi2201 demon- 
strate unambiguously that the ~70 meV nodal dispersion 
kink is due to electron coupling to the multiple phonons. 
In particular, we present a conclusive evidence that the 
electron-phonon coupling is strong in cuprates even in 
the heavily over-doped regime. 
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